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We map the conformation space of a simple lattice polymer chain to a network, where (i) the 
vertices of the network have a one-to-one correspondence to the conformations of the chain, and (ii) 
a link between two vertices indicates the possibility of switching from one conformation to the other 
by a single Monte Carlo move of the chain. We find that the geometric properties of this network are 
similar to those of small-world networks, namely, the diameter of conformation space increases, for 
large networks, as the logarithm of the number of conformations, while locally the network appears 
to have low dimensionality. 



The physical properties of polymers are the focus of a 
lot of attention due to (i) their role in many new advanced 
materials with important technologic applications [Q, 
and (ii) their role in biological processes [g| . An especially 
important unsolved problem in polymer studies is protein 
folding. Solving this problem has recently become more 
urgent as several studies suggest that a number of human 
diseases, including Parkinson's, Alzheimer's and British 
dementia, may be due to the aggregation of misfolded 
proteins 

The kinetics of protein folding are controlled by the 
structure of the free-energy landscape Theoretical 
calculations predict that in some cases the barriers in the 
free-energy landscape are quite small 1^. Hence, the dif- 
fusion of a protein's conformation on its energy landscape 
may be determined mainly by the structure and connec- 
tivity of conformation space. For this reason, much work 
has been done on the structure of conformation space and 
many models have been proposed, including tree struc- 
tures, random networks, and ultrametric spaces 0. In 
spite of all this work, very little is known about the con- 
formation space of polymers in general and proteins in 
particular. 

Here, we present evidence supporting the hypothesis 
that the conformation space of a lattice homopolymer 
chain may be a small-world network. These networks 
— which appear as the result of randomly replacing a 
fraction p of the links of a d-dimensional lattice with new 
random links — interpolate between the two limiting 
cases of a regular lattice [p = 0) and a random graph 
(j) = 1). The small- world regime is characterized by the 
properties (i) that a local neighborhood is preserved — as 
for regular lattices § — and (ii) that the average shortest 
distance between two vertices increases logarithmically 
with the number of vertices n of the network — as for 
random graphs [||. 

To gain some insight on the conformation space of real 
polymers, we consider here the conformation space of a 
2-dimensional lattice polymer chain. To study the ge- 
ometrical properties of conformation space, we map it 
onto a network . We first enumerate all allowed con- 
formations of the chain. We then identify (i) each confor- 
mation of the chain with a vertex of the network, and (ii) 



the possibility of changing from one conformation to an- 
other, through a single Monte Carlo move of the chain, 
with the existence of a link between the corresponding 
vertices (Fig. |l|). 

Our study relies on two important simplifications. The 
first simplification is to use a lattice model. On this re- 
gard, note that due to the limited number of allowed equi- 
librium angles between two monomers, the modeling of a 
polymer by a linear chain on a lattice may be experimen- 
tally justified [0 . In fact, even for the study of proteins, 
lattice models led to important insights (l^Jll|] . The sec- 
ond simplification is to neglect interactions. On this re- 
gard, note that here we are only interested in determining 
the geometric and structural properties of conformation 
space, and that interactions affect the rate of transfer 
between the allowed conformations |^,|l^,|l2|. Moreover, 
our results will also apply to any polymer chain in the 
limit of very high temperature for which the monomer 
interactions are mostly due to steric effects [p^ . 

Wc first calculate the geometric properties of confor- 
mation space for the simplified lattice polymer model and 
compare these properties with the predictions of several 
models for conformation space. This approach is differ- 
ent from the ones considered so far in the literature. We 
do not start by postulating a particular type of structure 
with some theoretically-desirable properties but, instead, 
try to generate the full network describing conformation 
space and compare its properties with different models. 

As a first test, we study the dependence of the aver- 
age shortest distance d between any two vertices in the 
network on the number n of vertices. Note (i) that the 
size of the network equals the number of allowed confor- 
mations of the chain, and (ii) that the distance between 
two vertices equals the minimum number of elementary 
moves of the chain necessary to switch between the two 
corresponding conformations. For a tree structure or a 
random network we expect a logarithmic increase of i 
with n, while for a c?-dimensional lattice we expect an al- 
gebraic increase: I ~ n^/'^ . In contrast, for a small-world 
network, I follows the scaling law jlij 

^(n,p)^(n*)i/''/(nK), (1) 

where the scaling function f{u) has the limits f{u) ^ 
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yi/d fgj. y_ <^ X ajj^cj ~ Inu for u ^ 1; n* ~ is a 
crossover size that separates the large- and small-world 
regimes, and p is the fraction of "rewired" links [Q. Fig- 
ure ^ displays our results for the polymer chain, which 
suggest that the conformation space can be described by 
a small- world network with p w 10^"^. Figure || clearly 
rules out the possibility that, for n 3> 1, the conformation 
space is a low dimensionality lattice. 

As a second test, we study the local structure of con- 
formation space and compare it with that of a random 
network To this end, we calculate the clustering co- 
efficient C, which is defined as the average ratio of the 
number of existing links between neighbors of a vertex 
and the maximum number of possible links. For a ran- 
dom network wc expect C ~ z/n, where z is the 
average connectivity of the network. In contrast, small- 
world networks have values of C of the same order of 
magnitude as those of regular lattices, because only a 
small percentage of links are different from those in the 
lattice In Fig. ^(a), we compare the values of C ob- 
tained for the networks with the values of C for random 
networks with the same size and connectivity. Clearly, 
the measured clustering coefficients are much larger than 
the expected values for random networks, ruling out a 
purely random structure for conformation space. 

As a third test, we calculate the number of elemen- 
tary loops in conformation space — usually referred to 
as the cyclomatic number |15| — and compare the re- 
sults with those for a tree structure. For a tree structure 
the cyclomatic number is identically zero while for all 
other networks it increases linearly with n. Figure ^(b) 
shows that the cyclomatic number for the polymer con- 
formation space clearly increases with n, consistent with 
a small-world network but ruling out a tree structure for 
conformation space. 

As a final test, we calculate the percentage of 
triplets {A, B, C} of vertices in conformation space whose 
distances obey an ultrametric relation dAC 5: 

max{dAB, dsc)- Figure ^(c) shows the percentage of ul- 
trametric triplets for conformation space. The percent- 
age of ultrametric triplets is significantly smaller than 
100%, ruling out a purely ultrametric structure for con- 
formation space. Moreover, the measured percentage of 
ultrametric triplets also rules out the random network 
and the tree structure as descriptions of conformation 
space. 

In summary, the regular lattice is rejected by the first 
test, the tree structure is rejected by the first, third, and 
fourth tests, and the random network is rejected by the 
first, second and fourth tests. Hence, we conclude that 
the geometrical properties of the conformation space of a 
lattice polymer chain are consistent with those of small- 
world networks but not with the geometric properties of 
the other geometries discussed in the literature. 

Next, we address the implications of our finding that 
the conformation space of a lattice polymer chain may be 
a small- world network. A central problem in protein fold- 
ing is the characterization of relaxation processes. The 



kinetics of a protein's conformations can be mapped to 
a random walk on conformation space. Naturally, the 
model studied here is too simplistic to enable us to un- 
derstand protein folding, however, the study of diffusion 
on conformation spaces with small- world geometries may 
still develop our understanding of relaxation in protein 
folding. Hence, we consider next the problem of diffusion 
on a small-world network. The usual way to solve a diffu- 
sion problem is to obtain the density of states p(A) of the 
transition matrix of the system |16|. For a small- world 
network, one has 16 17|] 



p{\) ^ A'^/^-i exp(-p/VA) 



(2) 



where the A's are the eigenvalues of the transition matrix 
for the network. To quantify the relaxation properties, 
we calculate the probability of return to the origin, usu- 
ally considered in the study of random walks |18|| , which 
is the Laplace transform of (0), 



dX p{\) exp(— At) , 



(3) 



where the lower limit in the integral is a cut-off that takes 
into account the finite size of the network. The return 
probability has the scaling behavior 



CXp(-(p2t)l/3) 

exp(-i/^^) 



t < ti 

ii < i < t2 , (4) 

t > t2 



where Po(oo) — l/n and the two crossover times scales 
can be written, in the small-world regime, as 



t2 



(n*)2(lnn)3 



(5) 



where we used the result £ ~ n* \nn |1J] . The results 
(^-|^) are consistent with our numerical simulations both 
for small-world networks and for the conformation space 
of the lattice polymer chain. 

The description of the conformation space as a small- 
world network may be useful for other complex disor- 
dered systems such as spin-glasses Specifically, 
small-world networks combine the features of apparent 
infinite dimensionality and low connectivity thought to 
be important for glassy relaxation ||l^,^. It has been 
hypothesized, that in order to reproduce stretched expo- 
nential decays the space of accessible conforma- 
tions must have a tree-like structure We have 
shown that for conformation spaces with a small-world 
structure stretched exponential relaxation can appear as 
an intrinsic feature of the space's geometry. 
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FIG. 1. The conformations of a lattice polymer chain and 
the mapping of the conformation space to a network, (a) 
The five panels labeled ci to cs show 5 different conforma- 
tions of a 15-mer constrained to have the end monomers fixed 
a distance A apart. We use circles with a light shade to rep- 
resent the fixed monomers, and dark circles to represent the 
moving monomers. For numerical convenience, we constrain 
the chain to occupy the half plane bounded by the thick black 
line. Note that the same results are obtained for other bound- 
ary conditions. The conformation of the chain evolves by the 
usual Monte Carlo elementary moves: the corner flip and the 
"crankshaft". For example, the chain can switch from con- 
formations ci to C2, C2 to C3, and C3 to C4 by single corner 
flips. In panels C2 to C5, we use dashed lines and dashed cir- 
cles to represent the position of the monomer moved from the 
previous conformation. We use conformations C4 and C5 to 
illustrate the "crankshaft" move, which involves the simul- 
taneous movement of two monomers, (b) Mapping of the 
conformation space of a chain to a network. We first allo- 
cate a vertex of the network to each allowed conformation of 
the chain. We then create a link between two vertices if the 
two corresponding conformations differ by a single elementary 
Monte Carlo move. For example, since we can switch between 
conformations ci and C2 by a single move, we put a link be- 
tween the vertices ci and C2. On the other hand, because C5 
cannot be reached from ci through a single elementary move, 
we do not place a link between the corresponding vertices. 
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FIG. 2. Average shortest distance between chain con- 
formations. We generate all the possible conformations of 
the polymer chain defined in Fig. |l| for A = 1,...,6 and 
for chains with m = 6, . . . , 15 monomers. For each pair of 
values (A,m), we identify in the conformation space all the 
different connected networks and calculate their size n. Here, 
we calculate the average shortest distance between every pair 
of vertices using to the breadth-first search algorithm, (a) 
Loglog plot of the average distance, rescaled according to the 
value of average connectivity [^. For n < 100, the average 
distance £ increases as n^^'^, as it would for a 2-dimensional 
lattice, (b) Same data but in log-linear plot. For n > 1000, 
we observe ^ Inn, as for a random network, and as pre- 
dicted by Eq. (jl|). Our results are consistent with the case 
of a small- world network with p ~ 10~^, clearly ruling out 
a low dimensional lattice as a model of conformation space. 
Note that for this and all other figures, the results for the 
small-world network have been average over 100 realizations 
of the network while each symbol "o" corresponds to a single 
conformation network. Hence, there is far less noise for the 
small- world networks data. 
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FIG. 3. Structure of conformation space. (a) 

Normalized clustering coefficients for the polymer chain, 
C(n, 2)/Crandom(n, z) ~ nC{n,z)/z, where z is the average 
connectivity. For a small-world network, the clustering co- 
efficient will be approximately a constant (for fixed p), and 
the normalized clustering coefficient increases linearly with 
n. Our results clearly rules out a purely random structure 
for conformation space since the normalized clustering coef- 
ficient for the polymer chain conformation space is orders of 
magnitude larger than the value for a random network. Note 
that the apparently logarithmic discrepancy between the two 
curves in the figure may be due to the logarithmic increase 
of z with n, which is shown in the inset, (b) Cyclomatic 
number. The cyclomatic number for a tree structure is iden- 
tically zero because there are no loops. For a generic network 
[ p^ , the cyclomatic number increases as nz. Hence, we plot 
the cyclomatic number for the polymer conformation space 
normalized by the connectivity z. As expected, we observe 
a linear increase with size n. (c) Percentage of ultrametric 
triplets. We calculate the number of triplets obeying the rela- 
tion dAC < max(dAB, dsc). For a purely ultrametric space, 
all triplets have distances obeying this relation. For 2-di- 
mensional square lattices and for small-world networks only 
slightly more than 2/3 of the triplets — which is the lower 
bound — obey the ultrametric relation. In contrast, for ran- 
dom networks and tree structures, we find a higher percent- 
age of triplets to be ultrametric. It is apparent that values 
for the random network and tree structure are not consistent 
with the measured values for the polymer conformation space. 
Note that each data point for the polymer chain corresponds 
to a single conformation network while the results for the 
small- world network represent an average over 100 networks. 



